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Abstract. The first optical proposal for the realization of the two-bit version of the 
Deutsch-Jozsa algorithm [D. Deutsch and R. Jozsa, Proc. R. Soc. London A 493, 
553 (1992)] is presented. The proposal uses Stark shifts in an ensemble of atoms 
and degenerate sources of photons. The photons interact dispersively with an atomic 
ensemble, leading to an effective Hamiltonian in atom-field basis, which is useful for 
performing the required two-qubit operations. Combining these with a set of one- 
qubit operations, the algorithm can be implemented. A discussion of the experimental 
feasibility of the proposal is given. 



PACS numbers: 03.67.Lx, 42.50.-p 



| On leave of absence from Physical Research Laboratory, Navrangpura, Ahmedabad - 380 009, India 



Two-bit Deutsch-Jozsa algorithm using an atomic ensemble 



2 



It is known that the quantum computers are able to perform certain specific jobs 
like factorization, searching a database etc. much faster than its classical counterparts. 
Several quantum algorithms have been developed to demonstrate the power of quantum 
computers. For example, Grover has proposed an algorithm to search a quantum state 
from an unsorted database of 2 n states pQ. This requires 0(2 n / 2 ) repetitions of certain 
unitary operations, demonstrating quadratic speed-up than the classical computers. 
Deutsch-Jozsa algorithm \2 i can be used to identify a certain function of n binary 
variables as constant or balanced by a single enquiry, whereas classically it would require 
up to 2 n ~ l + 1 enquiries. The key of this enormous speed of quantum computers lies 
in the so-called quantum parallelism, which enables one to apply the same unitary 
operation simultaneously on a number of basis states. 

There have been several proposals and experiments on implementation of these 
algorithms. For example, Grover's algorithm has been implemented using nuclear 
magnetic resonance (NMR) in bulk systems [3]. NMR jUH], ion trap .6], and linear 
optical [Zj techniques have been employed to implement the Deutsch-Jozsa algorithm. 
Recently, an optical scheme to realize the one-bit version of the Deutsch-Jozsa algorithm 
was proposed [Hj. To date, most of the works on Deutsch-Jozsa algorithm is for the case 
of a single qubit and we need to realize this algorithm for larger number of qubits, so 
that the full power of this algorithm can be realized. In this paper, we show how this 
can be achieved using Stark shifts in an ensemble of atoms. In this model, two freely 
propagating photonic qubits encoded in their polarization modes interact dispersively 
with an atomic ensemble. This leads to an effective Hamiltonian which along with 
several one-qubit operations enables us to implement the Deutsch-Jozsa algorithm for 
two qubits. We further discuss the experimental feasibility of our scheme using current 
technology. 

We start with the main features of the algorithm. This algorithm helps to identify 
whether a bivalued function fi(x) (which can take only the values or 1) of a variable x 
is constant or balanced. Here x is the decimal equivalent of the bit string x 2 , x n ) 
for n bits, i.e., x = Yh=i Xi2 n ~ l . Thus, x can take any of the 2 n values between and 
2 n — 1, for different binary combinations. If for half the values of x, the function takes the 
value and for the other half, it takes the value 1, then the function is called balanced. 
The function will be called constant if it assumes either the value or 1 for all values 
of x. In classical sense, it is easy to verify that at most 2 n ~ l + 1 evaluations of the 
function f%{x) for different values of x are necessary to determine whether the function 
is constant or balanced [see Table ^ . The Deutsch-Jozsa algorithm, on the contrary, 
requires only a single "evaluation" of the function for the same. To demonstrate such 
an immense power of this algorithm, we show in Fig. 1 the basic circuit 9J to perform 
the n-bit version of the Deutsch-Jozsa algorithm. It can be shown that if all the qubits 
of the system A are in |0) state, then the function fi(x) is a constant function. But if 
at least one of these n qubits is in the state |1), then the function can be identified as 
a balanced function. In what follows, we will demonstrate the two-bit version of the 
Deutsch-Jozsa algorithm in an ensemble of atoms interacting with a quantized field. 
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Table 1. Different possible functions showing their characterizations for n = 2. Here 
x is the decimal equivalent of the bit string (xi,x 2 ), e.g., x = 2 for two-bit input (1,0). 
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Figure 1. Basic circuit to perform the Deutsch-Jozsa algorithm for n qubits. All the 
n qubits in the system A are initialized to the state |0), while the qubit in the ancillary 
system B is initialized to the state |1). h A 1 h B - Hadamard transformations of the 
qubits of the systems A and B, Uf i - unitary operation operating on the system A+B, 
IV'o) etc- - states of the system A+B, x € {0, . . . , 2™ — 1}, y <G {0, 1}. The detector 
measures states of each qubit of the system A in (|0), |1)) basis. 
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Figure 2. A possible experimental set-up for the two-bit Deutsch-Jozsa algorithm. 
DOPA - degenerate optical parametric amplifier, PBS - polarization beam splitter, Ml 
etc. - lossless mirrors, HD1, HD2, VD1, VD2 - single-photon detectors for horizontal 
(H) and vertical (V) polarizations, HWP - half-wave plate, hi - sequence of quarter- 
wave and half- wave plates to perform the operation hi for photons, WP - sequence of 
quarter- wave and half- wave plates to perform the operations h! and h" for photons. 
The Hadamard rotations for the atomic ensemble arc not shown. 



Two-bit Deutsch-Jozsa algorithm using an atomic ensemble 4 



le> 




Figure 3. Relevant level configuration for implementing the Deutsch-Jozsa algorithm. 
The excited levels |e) and |e') are coupled to the ground levels \g) and \g'}, respectively, 
by the orthogonal polarization modes e+ and e_ of the photons, with equal detuning 
-A. 



In order to implement the Deutsch-Jozsa algorithm for two qubits, we need two 
distinguishable photonic qubits. Each photonic qubit is of course formed from two states 
of polarization. We would thus use two single photons propagating in opposite directions 
so as to keep the qubits distinguishable. Note that single photons have been used 
extensively as qubits for quantum computation based on linear optics [HI HI] . Further, 
the implementation of the algorithm requires a number of Hadamard transformations of 
photons which can be done by using quarter-wave and half-wave plates. We show in Fig. 
El a possible experimental set-up for the interaction of distinguishable photonic qubits. 
We use a degenerate optical parametric amplifier to produce two photons simultaneously 
|12| EJ- These photons interact with an ensemble of atoms. The experimental 
scheme shows how to make the photons distinct. The Deutsch-Jozsa algorithm requires 
that the polarization of two photons are measured simultaneously. For this purpose, we 
propose use of single-photon detectors [THj. The photon "1" ("2") can be detected in 
either of the detectors HD1 (HD2) and VD1 (VD2), which detects it in either horizontal 
or vertical polarization. The polarization beam splitter before the detectors preselects 
the polarization. Clearly, if the two detectors (HD1 or VD1 and HD2 or VD2) click 
simultaneously, the polarization of photons can be measured and thus the experiment 
leads to a conclusion for the algorithm. 

We consider the interaction of two distinct photons with an ensemble of N atoms 
where only four levels are relevant as shown in Fig. El Note that this kind of 
configuration can be found in the clock transitions in 133 Cs atoms and has been studied 
very extensively in ^01E|- Each photon has two circular polarization modes e±. We 
designate these modes as two orthogonal states 1 1 , 0} ^ and |0, 1)^, for the kth photon, 
where |1,0) (|0,1)) refers to a single photon in e + (e_) mode and no photon in e_ 
(e + ) mode. The photons in the states |1,0) and |0, 1) interact respectively with the 
l e ) ^ l#) an d <-» \g') transitions of each atom. We assume that the common 
frequency difference A between the photon polarization mode and atomic transition 
is much larger than the atom-photon coupling constant g, which is the same for both 
the modes. In this regime, the ground states \g) and \g') of each atom get Stark- 
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shifted. Using second-order perturbation theory, the respective shifts can be calculated 
as h\g\ 2 n + /A and h\g\ 2 n^/A, where n± are the total number of photons in e± modes. 
As the excited states |e) and \e') of the atoms remain unpopulated by interaction with 
the photons in large detuning regime, the effective Hamiltonian can be written as 

N 2 

#eff = h\J2 £ [\9j)(9j\\h0)k(l,0\ + W 3 )(9'M 1)*<0, 1|] , A = \g\ 2 /A .(1) 
j=i k=i 

We rewrite the above Hamiltonian in terms of the horizontal (|0)fc)and vertical (|l)fc) 
polarization states of each kth photon as 

i N 2 

H eS = 5 ftA£E[l&>fol(l >*<°l - |1>*<1| +<|0>*<1| +<|1>*(0|) 
Z j=l k=l 

+ |^)(^|(|0) fc (0|-|l) fc (l|-z|0) fc (l|-z|l) fc (0|)] , (2) 

where |0) fe = ^(|l,0) fc + |0, l) k ) and \l) k = ^(|l,0) fe - |0, l) k ) are identified as two 
basis states of the kth photonic qubit. As we deal with only two states of each atom, 
the ensemble of N atoms can have 2 N states. We use two of them as the states of the 
atomic qubit, namely YijLi \g)j{= |0)a) and ITjLi W)j{= • Note that as the atomic 
qubit consists of only the ground states of the atoms, decoherence due to spontaneous 
emission does not affect the process. Further, as the photons in free space can have long 
decay time, they do not lead to any effective decoherence as well. 

Under the action of the Hamiltonian (J2J), the following evolution occurs: 

|0) A |0, 0) F — > - e-* xm \0) A (-t\0) + 2 , 
\0) A \1, 1) F — > e- iXNt \0) A (-z\0) + 2 , 
|0) A |0, 1) F — - ze- aAf *|0) A (-*|0) + \l)r 2 , 

\0)a\1,0)f— > -ze~ aiV *|0) A H|0) + |l)f 2 , (3) 

where the state |0, l)p represents the photon 1 in the horizontal polarization state \0)p 
and the photon 2 in the vertical polarization state \l)p and so on. Similar evolution 
occurs when the atomic qubit is initially in the state |1)a- Note that the Hamiltonian 
(J2J) changes the polarization of each photon from linear to circular polarization, i.e., the 
atomic ensemble acts as a polarizer for each photon. We next show how this change in 
polarization helps us to implement the Deutsch-Jozsa algorithm for two qubits. Now 
onwards, we choose the interaction time T of the atomic ensemble with the photons 
such that XNT = tt/2. As shown previously jH], this interaction time is achievable 
with the available technology ^01 E|- For example, in an ensemble of 133 Cs atoms in 
a cell of dimensions 10x10x200 /im 3 at room temperature, the coupling constant g of 
each photon with the ensemble becomes 2.91 x 10 8 s" 1 and time T of interaction of the 
photon with the ensemble becomes 6.67 x 10" 13 s. Then for an atomic density of the 
order of 5 x 10 12 cm -3 (i.e., for iV = 10 5 ), the required detuning of each photon can 
be calculated from the above relation as A = 3.59 GHz m 12.35g. On the other hand, 
for Di transitions (S 2 ^!^ *-* 5 2 Pi/2, transition wavelength ~ 795 nm) in an ensemble 
of 87 Rb atoms in a magneto-optical trap (with diameter ~ 0.5 mm) ^H], we calculate 
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g ~ 3.53 x 10 6 s- 1 and T ~ 1.67 x 1(T 12 s. Thus for iV ~ 2.5 x 10 6 in the trap, we find 
A ~ 9g Clearly, the condition A ^> g is well satisfied in both cases. 

We now introduce the following four inequivalent Hadamard operations for each 
photon in (\0)k, |l)fc) basis and for each atom in the ensemble in (\g)j, \g')j) basis: 




I) 

1 -i 
-i 1 



(4) 



These operations can be implemented by using a resonant microwave field coupling the 
levels \g)j and \g')j of the jth atom. The corresponding interaction Hamiltonian is given 
by 



m = -nn 




g'h{g\ + a.c] , (5) 

where Q is the coupling constant of the microwave field with each atom in the ensemble 
and (f) is the phase of the field with respect to the matrix element between the levels 
\g)j and \g')j- Under the action of this Hamiltonian, the rotations (J3J) for each atom 
can be obtained by choosing Qt = tt/4 and = — ir/2, 0, ir/2, ir, respectively. Note 
that we do not provide the above rotations for the collective atomic states \0)a and 
\1)a- Alternatively, one could use two classical fields in Raman resonance to provide 
the equivalent rotation for each atom in the ensemble. For the photonic qubits, it is 
easy to verify that the operations (J3J) refer to rotation of linear polarizations to linear 
or circular polarizations, which can be performed using different sequences of quarter- 
wave and half-wave plates. Note that the matrix forms for transformations through 
quarter-wave plates and half-wave plates are given by 

cos(20) — i sin(20) \ 
sin(20) — cos(20) — i J ' 

( cos(20) sin(20) \ 

sin (20) -cos(20) ) ' W 

in (|0)fc, |l)fc) basis, where <fi is the angle of alignment of the wave-plate with the 
axis perpendicular with its plane. The single-qubit operations (j3J) for photons can be 
obtained using a series of these wave-plates of different orientations ["the SU(2) gadget"] 
as follows [1111201: 

h% = Q-K/iQ-n/iH-^/^ ; = Q-tt/4 ■ (7) 

Note that the half-wave plate rotates the polarization state from |0)fc to and vice 
versa for cf) = 7r/4. 

Further, as required by the Deutsch-Jozsa algorithm, two photonic qubits should be 
initially in the same polarization state |0)fc. This can be done by using a half-wave plate 
[see Fig. |2|. The degenerate parametric amplifier produces two photons in horizontal 
and vertical polarization states. One of them ("2") is in vertical polarization |1) and 
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thus is reflected by a polarization beam splitter. This photon is then sent through a 
half-wave plate (HWP in Fig- EJ) that rotates its polarization from |1) to |0), while the 
other photon ("1"), which is already in horizontal polarization state |0), passes through 
the polarization beam splitter. Thus both the photons are initialized to the state |0)fc. 

We now discuss how one can implement the Deutsch-Jozsa algorithm for two qubits 
using the evolution © and (jlj). We assume that two photons (the atomic ensemble) serve 
the purpose of the system A (B) in Fig. [TJ We start with the state j^o) f+a = |0, 0)^11)^ 
of the photons+ensemble system. Applying the Hadamard rotation h\ on each photon 
and the atomic qubit, we obtain 



F+A 



V2 



(|0> + |1» 



®2 TV 



J F 



\g')) ■ 



(8) 



Next the Uf t operations are applied. We show that in the present case, Uf t = h^Ueghf, 
where U e fi = exp[— iH ef[ t\ is the unitary operator corresponding to the Hamiltonian (J2J). 
Here the operations h eq correspond to different functions fi(x). 

In course of the operator sequence U^, the operation first prepares the atoms 
in (0)^4 state: 



F+A 



75(10) + ID) 



|o). 



J F 



Then the operation U e g for an interaction time T defined by XNT 
following [see Eqs. (J3J)]: 

52 



WD f+a = 



V2 



(-<|0> + |1» 



10} 



.4 



(9) 

7i / 2 yields the 
(10) 



^3,4 



We next identify the operations h eq for different possible balanced functions fi(x) [see 
Table as follows: 

-h" h' h cq = ti h" h eq = h" h" (M) 

- ilF 1 rl F 2 i "5,6 — a F 1 n F 2 > "7,8 — n F 1 a F 2 ■ \ Ll ) 

Here h! = hihjiz and h" = h\h<ihj, are composite sequences of different Hadamard 
rotations defined in Eq. (j3J) and yield the following: 

fc'|0) = e i7r/4 |0> , ti\\) = e" i7r/4 |l) ; 

3 i7r/4 |l) . (12) 



h"\0) = e" i7r/4 |0> , h"\l) 



Note that h! F . refers to h! operations on ith photon. The operators h eq lead to the 
following results, when applied on the state \ip") in Eq. (flDjl: 



±i(|0,0) + |0,l)-|l,0) 

±^(\o)-\i)) Fl ^=(\o) + \i)) F2 \o) A 



|i,i))f|o), 



V2' 



(13) 



l^2> 



/5,/6 



±i(|0,0)-|0,l) + |l,0) 



|i,i))f|o), 
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lv>2>/ 7 , /8 = w> 




1))f 2 |0) a 



(14) 



±^(|0,0)-|0,1)-|1,0) + 

± 4(l°> - |i>)^(|o> - 



1))f 2 |0) a . 



1,1))f|0)a 



(15) 



The final step of the algorithm is to apply the Hadamard rotation /if on 
each photon, which leads to the following outcomes: \ip 3 ) f 3 j 4 = |0, 1)f\0)a, \i J 3)f 5 ,h = 
\1,0)f\0)a, and \ip 3 )f 7 j a = |0, 0)f\0)a- Thus, after the final Hadamard rotations, either 
at least one of the photons remains in \0)f state, if the functions were balanced. 

On the other hand, in case of constant function f\, the Ufa operation is equivalent 
to the identity operation, while for the function f 2 , the Uf 2 operation is equivalent to the 
NOT operation on the atomic qubit. This can be implemented using microwave field 
for fit = tt/2 and = tt/2 [see the Hamiltonian (J5J) ]. Then, the output states 1^3) faj 2 of 
the photons+ensemble system become |1, l)i?ITj=i -^(\g)j ± W)j)- This means that if 
the functions were constant, both the photons can be detected in the state |1) f- Clearly, 
by measuring the polarization states of two photons at the end of the algorithm, one 
can characterize whether the functions were constant or balanced. 

To this end, we discuss possibility of decoherence in the entire process. As the 
photons interact with the ensemble dispersively, the excited states of the atoms get 
hardly populated. So the effect of spontaneous emission during the process is minimal. 
However, there could be collisional relaxation of the atomic ground states, which leads to 
loss of coherence in these states. But the time-scales of both these relaxation processes 
are of the order of 10~ 6 s, whereas, interaction time T of the photon with the ensemble 
is of the order of 10~ 13 s. Clearly, this is much less than the time at which the effects 
of different relaxation processes set in significantly. Thus, our model with photons and 
atomic ensemble is virtually decoherence-free. 

In conclusion, we demonstrated how the Deutsch-Jozsa algorithm can be 
implemented for two qubits using an atomic ensemble. One requires two distinguishable 
photons for this. We show how this can be achieved by using the output of a degenerate 
parametric amplifier. These photons dispersively interact with an ensemble of N atoms. 
The resulting Stark shifts of the atomic states lead to an effective Hamiltonian, which 
along with several Hadamard transformations, are used to implement the algorithm. 
We provide all the relevant operations to implement the algorithm and discussed the 
expected outcomes. 
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